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SUMMARY 


The atmosphere is acoustically dissipative and causes significant attenuation of 
propagating sound in addition to the attenuation caused by spherical spreading. Since 
the effect is frequency dependent, resulting in observed spectra that are not representa- 
tive of the source, it is important to account for it in much acoustical work. Commonly 
used acoustical data are the results of fractional- octave band analyses of spectra, 
whereas atmospheric attenuation is known fundamentally for discrete-frequency sounds. 
Correct methods of accounting for atmospheric attenuation in band data require consid- 
eration of the integrated effect across the bands for the specific distance involved. Such 
methods have not found their way into general use. Published analyses and computer 
programs have problems that inhibit implementation of correct methods by potential 
users. This report illustrates the underlying simplicity of an integral approach, docu- 
ments the basic analysis, and gives examples of errors incurred by not using correct 
methods. Computer programs are provided that are understandable, efficient, and 
simple to use. It is hoped that this will facilitate more widespread use of correct com- 
putational methods, especially where routine computer processing of data is employed. 


INTRODUCTION 

Spectrum analysis of air -propagated sound is fundamental to experimental acoustics. 
Most spectrum analysis data, particularly where they are associated with human an- 
noyance problems as in aircraft noise certification or in the general area of noise re- 
duction research, are comprised of sound pressure levels in 1/3 -octave frequency bands. 
These wide -band data are commonly used as an alternative to narrow -band or power 
spectral density data to allow frequency-dependent calculations to be made with a rea- 
sonably manageable am omit of data. 

It is well known that the atmosphere is acoustically dissipative and causes significant 
attenuation of propagating sound in addition to the attenuation caused by spherical 
spreading. The effect varies with frequency so that, as propagation distance increases, 
the observed spectrum differs more and more from that emitted by the source. Thus, 
a pervasive problem in noise analysis is to account for atmospheric attenuation - either 
to remove the effect in order to obtain so-called lossless data, which reveal the true 
character of a source, or to incorporate it into lossless data to make practical propaga- 
tion calculations. Often it is necessary to adjust measured data to standard conditions, 



since atmospheric attenuation is variable with ambient conditions. This requires two 
steps: removal of the effect due to the measurement conditions, and then reincorpora- 
tion of it for standard conditions. 

Atmospheric absorption is defined fundamentally for discrete-frequency sounds, and 
is usually expressed as an attenuation per unit distance. Analytically it can be specified 
as a continuous function of frequency and atmospheric variables. The nature of the 
phenomenon is well understood, and there is a body of literature on the subject. For 
this report, the physics of the phenomenon is not an issue. It is presumed, simply, that 
a suitable analytical representation of discrete-frequency atmospheric attenuation is 
available in the form of a continuous function of frequency. 

When the detailed spectrum of a sound is known, as in having the power spectral 
density function or the results of narrow-band analysis, accounting for the effect of 
atmospheric absorption is trivial in principle. The appropriate calculations are made 
either functionally or discretely over the frequency range. In practice this is rarely, if 
ever, done. 

For the more commonly used data obtained from wide-band analysis of spectra, 
however, the problem is less simple. Since a band sound pressure level is by defini- 
tion a measure of the integral of the power spectral density between the lower and upper 
band- edge frequencies, a single -valued measure of atmospheric absorption for a band 
can only be obtained by considering the integrated effect over the band. Occasionally 
this fact is not recognized. Most often, however, it is, but for the sake of expedience 
the rigorous calculation called for is deliberately avoided. Instead, the band attenua- 
tion is taken as that of some discrete frequency, the value of which is determined by 
some rule of thumb. 

The preeminent example of this approach is the widely used procedure specified in 
ARP-866A (ref. 1). Attenuation for a frequency band is normally obtained by using the 
discrete frequency attenuation at the band center frequency. For cases where the slope 
of the spectrum is suitably conditioned and at low frequency where the slope of the at- 
tenuation function is small, this is a reasonable approximation. Unfortunately, these 
conditions do not always prevail, especially at high frequency. This is partially ac- 
counted for in ARP-866A by the requirement that the lower band-edge frequency should 
be used instead of the center frequency when the latter exceeds a certain value. But the 
procedure still is approximate and requires suitably conditioned spectra if errors are 
to be small. The ARP-866A procedure also results in band atmospheric attenuations 
that are linearly proportional to distance. It will be shown later that this also is an ap- 
proximation. 

A correct computational method to obtain single-valued band attenuations requires 
consideration of integrals across the bands. However, such correct methods have not 
found their way into common use, and approximate methods prevail. This is true even 
where routine computer processing of data is employed, which seems to be unjustifiable. 
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Errors of varying magnitudes are introduced that depend on the shapes of the spectra 
involved, the absolute value of the frequencies, and the propagation distances. Ex- 
amples of these errors are given later. 

Recent publications (refs. 2 and 3) deal with the correct computation of band at- 
tenuations. However, in the author’s view, the analyses given therein are neither as 
succinct nor as understandable as they might be. Furthermore, the computer programs 
that result are lengthy and complicated, and the methods have other problems that are 
discussed later. These conditions would seem to keep potential users from implement- 
ing correct methods. 

The intent of the present work is to document the underlying analysis in a straight- 
forward way and to develop a method incorporated in computer programs that are simple, 
understandable, and easily employed. This is done in the hope that correct methods will 
find wider routine use, especially where computer processing of data is employed. 

The theoretical development given herein is rigorous and makes use of the relation- 
ships, which are a matter of definition, between the power spectral density of a sound 
signal and wide-band analysis representations of it. Since the resulting method is in- 
tended for computing band atmospheric attenuation when only band levels are known, the 
spectrum shape within each band must be inferred from the band data. This is the only 
substantive approximation required. However, rather than being a weakness, this ap- 
proximation is a key element that distinguishes the integral approach from the discrete- 
frequency method. Any alternative scheme for inferring the spectrum shape may be 
substituted at the discretion of the reader. The one given herein, which uses straight- 
line segments, has been used satisfactorily by the author since 1968. 

Computer programs, one for performing the band integrations and another for gen- 
erating discrete-frequency atmospheric attenuation, have been published previously 
(ref. 4). However, a detailed analysis of the procedure was not given, and, in fact, 
the band integration program contained some shortcomings. These included a lack of 
distinction between lossless or air -attenuated input data, a linear dependence of band 
attenuation on distance, and a minor frequency -stepping error across each band. These 
are all removed here. 

The discrete -frequency attenuation program of reference 4 was based on the graph- 
ical information given in reference 5. It should be noted that the specific discrete- 
frequency attenuation model employed is not critical to the development of the band in- 
tegration procedure, and any model may be substituted. In reference 2, an improved 
discrete -frequency attenuation model is presented. This model has been adopted for the 
present work on the basis of its merits, particularly the wide frequency range of its 
applicability. 

The band attenuation computational procedure developed in reference 2 also has 
some problems. The most significant is that the number of bands diminishes by two 
each time an atmospheric attenuation calculation is made. The first and last bands are 
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lost because of the computational technique employed. Thus, if air-attenuated data are 
to be converted to standard-day data through the intermediate step of obtaining lossless 
data, which is a common requirement, a total of four bands will be lost. 

The method developed herein eliminates this shortcoming. In addition, the present 
development is structured with a clear separation between the numerical integration 
procedure and the technique to infer the spectrum shape from the band data. This makes 
it easy for the reader to modify or replace the technique used to infer the spectrum 
shape. Furthermore, the formulations are developed to facilitate computer coding, and 
frequency incrementing is accomplished by methods that generate so-called preferred 
frequencies. 

Finally, the development given herein defines a band as that portion of a spectrum 
between certain frequency limits. This is equivalent to specifying measured data ob- 
tained with ideal filters. In reality, band data are obtained experimentally with band- 
pass filters. It has been shown in reference 2 that for extreme spectrum or attenua- 
tion gradients the filter skirt can have a significant effect on a spectrum analysis. Such 
extreme conditions are not usual, however, and to account for the effect of the filter 
skirt greatly complicates the development. The attendant discussion would need to go 
far beyond the scope of the present work. For the usual frequency range from 50 to 
20 000 hertz, and even to higher frequencies as will be seen, the filter skirt is not 
likely to present a problem. For further discussion of this point, and for quantitative 
estimates of the conditions under which filter characteristics can present problems, the 
reader is referred to reference 2. 


SYMBOLS 

Aj(r) band atmospheric attenuation, dB 

a(f) discrete -frequency atmospheric attenuation per unit distance, dB/m 

b fractional- octave band size designation denoting 1/b-octave bands 

C constant 

f frequency, Hz 

f band center frequency, Hz 

fp f 2 lower and upper band-edge frequencies, respectively, Hz 
g function 

i band index 

j subband index 
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I 

incremental propagation path index 

lossless fractional- octave band sound pressure level, dB 
air -attenuated fractional-octave band sound pressure level, dB 

to eq. (14), lossless subband sound pressure level; for eq. (21), lossless or 
air -attenuated subband sound pressure level, dB 

air -attenuated subband sound pressure level, dB 

subband sound pressure level relative to subband level at band center fre- 
quency, dB 

integer 

band number 

number of subbands in fractional- octave band 
subband rms acoustic pressure, Pa 
propagation path length, m 
sum 

power spectral density of lossless acoustic pressure, Pa /Hz 

2 

power spectral density of air-attenuated acoustic pressure, Pa /Hz 

ANALYSIS 

Let W(f, r) denote the power spectral density of the rms acoustic pressure of freely 
propagating sound at any distance r from a source in the absence of atmospheric at- 
tenuation. This represents the source (lossless) spectrum. If a(f) denotes the discrete- 
frequency atmospheric attenuation in units of decibels per unit distance, the air- 
attenuated spectrum W’(f, r) in a homogeneous atmosphere is expressed as 

-a(f)r 

W’(f, r) = W(f, r)10 10 (1) 

(It is understood that atmospheric attenuation is dependent on ambient temperature and 
humidity. These variables are omitted from the functional notation a(f) for simplicity. ) 

The sound pressure level L^r) of the portion of a spectrum in the i^ band be- 
tween the frequency limits of (fj). and (f£)j is expressed in terms of the power spectral 
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density as 


ft 2 >i 

L^r) = 10 lflgy (fl) W(f, r)df + C (2) 

where C is an arbitrary constant required by the definition of a level. * The air at- 
tenuation in decibels A.(r) for a spectrum band is the difference between the band level 
of the lossless spectrum at some location r and the corresponding band level of the 
air -attenuated spectrum: 


/*(f2). Wf2). 

A.(r) = 10 log 7 (fi) 1 W(f, r)df - 10 log 1 W’(f, r)df (3) 

If the air -attenuated spectrum in the second term on the right of equation (3) is re- 
placed by equation (1), equation (3) becomes 



-a(f)r 

W(f, r)10 10 


(4) 


Equation (4) expresses the band attenuation in decibels in terms of the lossless spec- 
trum. 

It is just as likely that the air -attenuated spectrum and not the source spectrum will 
be known. It will be necessary to compute the band attenuation from these data also. 

To develop a suitable expression for this case, equation (1) is rewritten as 

+a(f)r 

W(f, r) = W f (f, r)10 10 (5) 


•*As noted in the INTRODUCTION, this formulation is equivalent to specifying a 
measurement using an ideal filter. If a filter transmission characteristic t(f) is to be 
accounted for, where 0 < t(f) < 1, the function t(f) must be introduced as a factor into 
the integrand and the limits taken from 0 to °°. Subsequent development would proceed 
accordingly. 
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and substituted in the first term on the right of equation (3) to obtain 


/ (f 2>i +a(f)r 

W’(f,r)10 10 df - 10 log I W’(f, r)df (6) 

x‘l> 4 •'(fl)i 

Equation (6) is equivalent to equation (4) but is expressed in terms of the air -attenuated 
spectrum. 

Note that equations (4) and (6) both contain the propagation distance r within an in- 
tegrand from which it cannot be removed. This establishes a nonlinear dependence of 
band attenuation on distance, which discrete -frequency approximate methods inherently 
neglect. 

It is also appropriate to note here that the foregoing analysis can be simply ex- 
tended to an inhomogeneous atmosphere. The quantity a(f)r expresses the total 
discrete -frequency attenuation in decibels along the propagation path. To extend the 
results for an inhomogeneous atmosphere it is only necessary to define the total 
discrete-frequency attenuation along the path as 

X a ( f > dr 

where a(f) is an implicit function of r by virtue of its dependence on atmospheric vari- 
ables. Such an extension for making practical calculations is considered later. 

If the symbols L^(r) and L!(r) are used in place of the integral expressions for the 
band levels in equation (3), simple rearrangement results in 

L[(r) = L i (r) - A.(r) (7) 

and 

L^r) = L!(r) + A^r) (8) 

Lossless, freely propagating sound obeys the inverse square law so that, for two dis- 
tances r j and r ^ along the path of propagation, the lossless spectrum band levels are 
related by 



If r 2 replaces r in equation (7) and equation (9) replaces the first term on the right, 
equation (7) becomes 


L!(r 2 ) = L^rj) + 20 log^— ^ - A^rg) (10) 

Also, substituting r^ for r in equation (8) and replacing the first term on the right of 
equation (9) with the result yield 


L i (r 2 ) = L!(rj) + 20 log^ + A.frj) (11) 

Equations (10) and (11), in conjunction with equations (4) and (6), respectively, pro- 
vide the means to account for atmospheric attenuation in spectrum analysis band data. 
Equations (4) and (10) permit determination of the air -attenuated band levels at any dis- 
tance when the lossless levels at some distance are known. Equations (6) and (11) per- 
mit determination of the lossless band levels at any distance when the air -attenuated 
levels at some distance are known. 


NUMERICAL APPROXIMATION 

Equations (4) and (6) contain integral expressions involving power spectral density. 
It is obviously easy to obtain equivalent finite-sum expressions that can be evaluated as 
precisely as required. However, in most practical computational problems, only band 
sound pressure level data are available - not power spectral density, which the expres- 
sions apparently require. Fortunately, there is a natural connection between power 
spectral density and the spectral increments involved in finite sum expressions, which 
permits the power spectral density to be replaced by sound pressure levels in finite fre- 
quency subintervals within a band. These sound pressure levels can be obtained by in- 
ference from the band data, and essentially reflect the spectrum shape. In this section, 
numerical integration expressions for band attenuation and related equations for fre- 
quency incrementing are developed, and the scheme to infer the spectrum shape from 
the data is described. 
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Band Attenuation 


The integral expressions for band levels can be approximated directly by finite 
sums. Equation (2), for example, can be approximated by 

n 

L t (r) = 10 log ^ W(fj, r)(Af) j + C (12) 

j = l 

th 

where the i frequency band is divided into n contiguous intervals, each embracing fj 

and each of width (Af).. (Note that for notation simplicity the subscript i is omitted 

J 

from the frequency quantities on the right since their restriction to the i band is 
clear. ) The representation may be made as accurate as necessary by increasing the 
value of n and choosing fj and the intervals (Af)^ in such a way that W(fj, r) approxi- 
mates the mean value on each interval. 

By definition, the power spectral density is related to the rms acoustic pressure 
Pj(r) in a band embracing L by 


w(f j9 r) 


p ft r) 

< i( >i 


< A ,) r 0 


Thus, equation (12) may be written 


L.(r) = 10 log ^ pf(r) + C 

j = l 


(13) 


where the p.(r) (j = 1, 2, . . . , n) denote the rms acoustic pressures in sufficiently 
j th 

small intervals that collectively comprise the i band. 

From the foregoing and the fact that the rms sound pressure is related to the sound 

pressure level Lj(r) decibels (which is the customary form for acoustical data) by 


P] 2 (r) 


Pref 


10 


10 
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equations (4) and (6) can be written as 



Equations (14) and (15) express atmospheric attenuation, in decibels, for bands of 
frequency as functions of sound pressure levels in smaller intervals within the bands. 

To carry out computations it is only necessary to choose frequency intervals across the 
band and to determine the sound pressure levels in the intervals. 


Frequency Intervals 

For fractional- octave band analysis, the bandwidth increases in proportion to fre- 
quency and the data are conveniently exhibited on a logarithmic frequency scale. It fol- 
lows that, for summation purposes in equations (14) and (15), sound pressures are more 
appropriately specified at equal increments of the logarithm of frequency. 

Delineation of a fractional- octave band on a logarithmic frequency scale is shown 
schematically in figure 1. The band center frequency f is shown by the bold dashed 
line and the center frequencies of consecutive bands would appear at equal spacing. The 
frequencies of the lower and upper band-edges are shown by the bold solid lines. The 
center frequency of a band is the geometric mean of the lower and upper band-edge fre- 
quencies (i. e. , f = y/Tfo), and on a logarithmic scale it appears at the center of the 
band. 

The band is shown subdivided by light solid lines into n frequency intervals con- 
sistent with the requirements for the summations of equations (14) and (15). Herein- 
after these intervals are referred to as subbands. The center frequency of each sub- 
band is denoted by the short dashed lines. As will be seen, it is convenient to use an 


10 



odd number of subbands; therefore, n is defined according to 


n = 2m + 1 


( 16 ) 


where m is some integer. The subbands can now be indexed on j from -m to m, as 
shown in the figure, resulting in the subscript j=0 for the center subband. 

If the notation 1/b-octave is used to designate the size of fractional-octave bands, 
where b is some number (usually, but not necessarily, an integer), the center fre- 
quencies of successive bands are obtained from 


(f c )j = 10 10b (17) 

as i takes successive integer values. (This formula restricts frequencies to standard 

preferred values as specified in ref. 6.) Referring to figure 1, the reader can verify 

th th 

that the center frequency of the j n subband in the i band is given by 


( f j>i 


«c 



(18) 


Equation (18) permits regular indexing of frequency across the band for use with equa- 
tions (14) and (15) if the summation limits go from -m to m. 


Programming Considerations 

Equations (14) and (15) with indexes revised according to the discussion of the pre- 
ceding section each involve summations S having a functional form 


S = 


m 

2 g(j) 

j=-m 


This can be written equivalently as 


S = g(0) + 



(19) 


( 20 ) 
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Although the form of equation (20) is less simple than that of equation (19), equation (20) 
is more suited to programmed computation. For either form, an initial value for the 
accumulating sum must be set. However, equation (20) allows the initialization step to 
be used to compute g(0). The summation of equation (20) then proceeds from 1 to m, 
which requires only half the number of cycles than equation (19) requires. The sub- 
script j=0 denotes the band center frequency; hereinafter the 0 is replaced with c. 

If the foregoing considerations are applied to equations (14) and (15), the function 
of r notation is dropped for simplicity, and some rearrangement is made, these equa- 
tions can be combined into 



where L is now used to denote either lossless or air-attenuated subband levels. The 
plus signs are used with air -attenuated levels and the minus signs are used with loss- 
less levels. 

Note that a constant added to all levels in equation (21) has no effect on the result 
since it introduces a common factor into all terms of both numerator and denominator. 
This permits equation (21) to be further simplified by expressing all levels relative to 
the level at the band center frequency. Equation (21) now can be written in final working 
form as 



( 22 ) 
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where AL is the subband level relative to the subband level at the band center fre- 
quency. Values of frequency for use in equation (22) are obtained from equation (18). 

It was noted earlier that extension of the analysis to an inhomogeneous atmosphere 
involved only the total discrete-frequency attenuation along the propagation path. This 
is the quantity a(f)r, which has retained its identity throughout. To use equation (22) 
for the inhomogeneous case, it is only necessary to replace the total attenuation along 
the path with 

2 a k( f H Ar )k 

k 

where the subscript k denotes subdivision of the path into small increments Ar over 
which a(f) is constant, and for which 


2 (Ar) k = r 
k 


Such an extension is well suited to a layered atmosphere model. 


Subband Sound Pressure Levels 

It is not possible to obtain the rms pressure distribution across a frequency band 
from band data alone, since each band level inherently represents the integration of the 
distribution. However, the distribution can be approximated. 

It can be shown that a power spectral density function that varies as a constant 
power of frequency (i. e. , the power spectral density level exhibits a straight line re- 
lationship on a log frequency scale) results in fractional-octave band levels that likewise 
vary as a constant (but different) power of band center frequency, irrespective of the 
size of the bands. This establishes that, for fractional- octave band spectrum analyses 
presented graphically as band pressure level against logarithm of frequency, the slope 
of the line connecting the band levels at their respective center frequencies is inde- 
pendent of the size of the bands. The absolute levels obtained from analyses using dif- 
ferent bandwidths, however, differ by ten times the logarithm of the bandwidth ratio. 
These properties are illustrated in figure 2. Thus, subject to an assumption about the 
smoothness of the spectrum, it is possible to infer from a fractional- octave band spec- 
trum analysis the analysis in yet smaller fractional- octave bands. The number of sub- 
bands into which a band may be divided may be chosen without limit. 
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On the basis of the foregoing considerations, a scheme has been developed to ap- 
proximate the spectral distribution across fractional-octave bands so that band atmos- 
pheric attenuation can be computed with the equations developed. This is illustrated in 
figure 3. Conceptually, for a given spectrum analysis, the levels of the fractional- 
octave bands, plotted at the band center frequencies on a logarithmic frequency scale, 
are connected by straight-line segments. The distribution (but not the absolute values) 
of the subband levels across each band from edge to edge is thus established approxi- 
mately by a pair of straight-line half-segments whose slopes are determined by the 
data. For the first and last bands, where only one of the half-segments is defined by 
the data, the slope is assumed the same on both sides of the band. 

It was shown in the previous section that only the relative distribution of the sub- 
band levels across the band was significant. Equation (22) therefore was written in 
terms of levels relative to the level at the band center frequency. Definition of the rel- 
ative subband levels using straight-line half-segments is illustrated in figure 4. The 
fact that only the relative distribution of the subband levels is significant restricts at- 
tention to the level increment from band to band; this is loosely referred to as the slope 
of the data. 

Finally, it should be obvious that when discrete tones are present in a spectrum in 
small numbers to the extent that one, or at most a few, influences the level of any par- 
ticular band the procedure given here does not yield an accurate subband level distribu- 
tion. In such a situation, where a tone -dominated band stands out and where for lack 
of information about the tone frequencies the procedure of this section is used straight- 
forwardly, the computed band attenuation value will be strongly influenced by the dis- 
crete tone value at the band center frequency. The error introduced is potentially of no 
more seriousness than neglecting the integral approach altogether. On the other hand, 
if the discrete tone frequencies (and their relative levels, if there are more than one 
within a given band) are known, accurate calculations can be made. In this case, the 
frequencies used for the summation of equation (22) would be obtained not by equa- 
tion (18) but by the tone frequencies known to exist, and the summation would proceed 
only for as many terms as there are discrete tones, possibly only one. If this were 
done, care would also be needed to remove the effect of the tones before establishing the 
line segments, and ultimately the attenuation, for the bands adjacent to the tone- 
dominated band. 


COMPUTER PROGRAMS 

Equations (10), (11), (17), (18), (22), and a description of atmospheric absorption as 
a function of frequency constitute a rigorous numerical method to calculate atmospheric 
attenuation in spectrum bands. The accuracy is limited only by the number of terms in 
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the summations and by the scheme used to infer the subband levels from the band data. 
The scheme used herein was described in the preceding section and is illustrated in 
figure 4. 

A pair of Fortran computer pregrams that embodies this method is given in the ap- 
pendix. The first pregram, BASPAT, which obtains the subband levels and performs 
the numerical integration, has been designed for efficiency, computational power, and 
simplicity of use. Input data comprise sound pressure levels in contiguous fractional- 
octave bands at some specified distance. The size of the bands is selectable. 

If the input data are lossless, air -attenuated data will be generated. If the input 
data are air -attenuated, lossless data will be generated. The distance for which the 
output data are obtained may be different than the distance for the input data. To adjust 
measured air -attenuated data to standard-day air-attenuated data, two calls of the pro- 
gram are necessary. The first call will generate lossless data, which must then be 
used in a second call to obtain standard-day data. 

The values of all frequencies used in the program are obtained by methods that 
generate preferred numbers consistent with international standards. For further in- 
formation on these standards, the reader is referred to references 6 and 7. Briefly, 
the center frequency of any band can be expressed as 


f 


c 


10 


N h /10 


(23) 


where N b is defined as the band number. The band number, moreover, must conform 
to 


N b = - n (24) 

D b 

where b denotes the kind of band as defined in the Frequency Intervals section and n 
is some integer. (Note that the band number need not be an integer. ) 

In using the program, the user specifies the values of the band center frequencies 
corresponding to the data simply by specifying a single variable. This is the interna- 
tional band number of the first band, denoted IBNF in the program. The user must in- 
sure that IBNF is in conformance with equations (23) and (24). The program uses 
this value to compute the center frequency of the first band and all succeeding bands of 
the air attenuated data. 

An alternative version of BASPAT is also given in the appendix. This is BASLAT 
and is functionally the same as BASPAT except that it allows calculations for a layered 
atmosphere in accordance with previous discussion. It is provided as a convenience to 
the reader, and all discussions herein with regard to BASPAT apply to it as well. 
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The program BASPAT calls a subroutine to obtain discrete-frequency atmospheric 
attenuation. The second program, ATMAT, is provided for such calls. The program 
ATMAT computes discrete-frequency atmospheric attenuation in accordance with the 
data presented in reference 2. Obviously, the reader can substitute his own version of 
this subroutine. 


SAMPLE CALCULATIONS 


In this section, results obtained by the method and computer programs given herein 
are obtained for comparison with results from the numerical integration method of ref- 
erence 2. In the subsequent discussion, examples of errors caused by discrete- 
frequency approximations of numerical integration methods are presented. 

To illustrate the use of the programs, and for convenient evaluation of the results, 
the sample problem of reference 2 is repeated here. The original data for starting 
the calculations appear in table I, column 3. The data are presumed to be lossless 
sound pressure levels 20 meters from the source. The first calculation gets the corre- 
sponding air -attenuated data at the same distance for an ambient temperature of 293 K 
and a relative humidity of 50 percent. To perform the calculation, the subroutine 
BASPAT is executed from a suitable calling program in which the input variables in the 
subroutine vector are set as follows: 


SLI array of original lossless 1/3-octave band sound pressure levels as shown in 
table I, column 3 

RI 20. 0, propagation distance from the source in meters associated with SLI values 

IBNF 36, which, in accordance with eqs. (23) and (24) for b = 3, specifies that the 

frequency associated with the first band level is 4000 Hz nominal (see table I, 
columns 1 and 2) 

NF 15, number of bands 

NB 3, designates 1/3 -octave bands 

T 293.0, ambient temperature, K 

RH 50. 0, relative humidity, percent 

IDK 1, declares input data are lossless, which will generate air -attenuated results 
RO 20. 0, distance from the source in meters for which the results are desired 

The results of this calculation by the program BASPAT are shown in table I, 
column 4. The corresponding results from reference 2 appear in column 5. As in ref- 
erence 2, these re; 'Us (column 4) are next used as data for the inverse calculation. 
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For this case, the variable IDK in the BASPAT calling vector is set to zero to declare 
that the input data are air attenuated. All other data remain the same. The results of 
the inverse calculation appear in column 6. The corresponding results from refer- 
ence 2, computed from column 5, appear in column 7. 


DISCUSSION 

Table I shows that the air -attenuated results from the two methods agree exactly 
at 6300 hertz. Below 6300 hertz, because of the reduced gradient of the discrete- 
frequency air -attenuation function with frequency, exact agreement prevails. Above 
6300 hertz, the results differ by a slightly increasing amount, never exceeding 0. 7 dec- 
ibel up to 80 000 hertz in spite of absolute adjustments to the original data of close to 
55 decibels. The differences are attributed to neglecting the filter characteristic in the 
present procedure. The filter characteristic is of significance where the gradient of 
the product of air absorption and power spectral density is high, which generally occurs 
at high frequency. The gradually increasing effect is due to the significantly increasing 
slope of the air attenuation. The filtered values, column 5, are higher because of the 
filter skirts, which pass energy from adjacent bands. Nevertheless, neglecting the 
skirt- results in fairly small errors, particularly in consideration of the extremely high 
frequencies of the example. The present method also results in no loss of data, which 
is not true of the method of reference 2. 

The inverse results obtained by the present method agree to one decimal place ex- 
cept at the highest frequency. A corresponding evaluation of the method of reference 2 
cannot be made because of the units precision reported for the lossless results (col- 
umn 7) and the loss of band data. 

With regard to obtaining lossless data from air -attenuated data accounting for the 
filter characteristic, it must be pointed out that reference 2 contains an implicit error. 
The development given therein postulates that the air -attenuated (received) spectrum is 
known. This spectrum corresponds to air -attenuated band data obtained from ideal fil- 
ters. However, the practical situation for which the equation is developed provides 
air -attenuated data obtained from real filters. Thus the resulting equation, though 
intrinsically correct, is not strictly applicable. Moreover, it is not possible to obtain 
ideal-filter data from real-filter data and the filter characteristic. A complete discus- 
sion of this problem is beyond the scope of this report. The fact that the inverse calcu- 
lation of reference 2 seems to be correct is attributable to the units precision of the re- 
sults, since the errors are expected to be less than 1/2 decibel for this example. 

Data are given next to illustrate the errors involved in not employing an integration 
procedure to account for air attenuation in band data. The errors shown in figure 5 are 
introduced by using discrete frequency attenuation at the band center frequency in place 
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of the attenuation obtained from the present numerical integration method. These 
examples illustrate the errors that will result in lossless data obtained from air- 
attenuated data for ambient conditions of 298. 15 K and 70 percent relative humidity. 

The dependence of the error on frequency and slope of the spectrum as designated by 
the data slope is shown. Note that these results are illustrative only and that errors of 
different magnitudes will occur for other ambient conditions and propagation distances. 
These data illustrate, however, that for short propagation distances, even to extremely 
high frequencies, the errors tend to be small. On the other hand, there are conditions, 
especially for longer distances, under which the errors can be significant. These ob- 
servations are consistent with the conclusions of reference 2. 

Figure 6 shows the errors that can result in lossless data by assuming that band 
attenuation is linearly proportional to distance. Here, too, the errors shown are in- 
tended to be illustrative, and they will differ for other conditions. 

The significance of these data is that under some conditions approximate methods 
produce errors that may exceed the precision required for various purposes. Obvious- 
ly, for suitably conditioned data these errors can be negligible. But to insure small 
errors requires that data be monitored. The potential for error of any significance can 
be fully eliminated along with the need for data monitoring by the routine use of numer- 
ical integration methods. 


CONCLUDING REMARKS 

The use of discrete-frequency atmospheric attenuation at some characteristic fre- 
quency in a band to represent the band attenuation is inherently erroneous. A method 
has been developed herein that properly computes fractional-octave band atmospheric 
attenuation. Ideal filters are assumed, and the total method is reduced basically to two 
short computer programs. Three logically separable elements are involved: a scheme 
to infer the spectrum shape from band data, a numerical integration procedure, and 
discrete-frequency atmospheric attenuation. The first two elements are incorporated in 
a short Fortran program of which two versions are given - program BASPAT, which is 
restricted to a homogeneous atmosphere, and program BASLAT, which is essentially 
the same but permits calculations for a variable atmosphere. This numerical integra- 
tion program, in both versions, has been written to be understandable and computa- 
tionally efficient. The discrete-frequency attenuation computational method is essen- 
tially taken from the literature and is given as the short program ATMAT. This 
program may be replaced at the discretion of the reader. 

It is shown by example that the effect of the filter is small and that the results com- 
puted agree with another published numerical integration method. Examples are also 
given of errors that occur by not using such methods. 
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The present method overcomes the shortcomings of other published methods, is in- 
herently simple, and is incorporated in short uncomplicated computer pregrams. In 
view of the correctness of the numerical integration procedure and the inherent sim- 
plicity of the present method and its implementation, this method is offered as a re- 
placement for discrete-frequency approximate methods where routine computer proc- 
essing of data is employed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 1, 1978, 

505-04. 
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APPENDIX - COMPUTER PROGRAMS 


Given herein are the Fortran computer programs that constitute the complete 
method to compute the effect of atmospheric absorption in sound pressure level data 
comprising the results of fractional- octave band analyses of propagated sound pressure 
spectra. The program listings contain descriptive information on their functions and 
use and on the variables employed. 


NUMERICAL INTEGRATION 


Subroutine BAS PAT 


Subroutine BASPAT (BAnd SPectrum ATtenuation) uses the scheme previously de- 
scribed to infer the subband sound pressure levels from the input data and performs 
the calculations specified by equations (10), (11), (17), (18), and (22). 

The value of M is coded as 3, which by equation (16), specifies that each band in- 
tegration is approximated by a sum of seven terms. This has been found by the author 
to be quite satisfactory in general, but it may be easily changed at the discretion of the 
reader by changing the value of M. (If M^is made greater than 5, it will be necessary 
to increase the dimension of CL and ) d 

The variables SLl and SL2 are the values of the subband levels used in the sum- 
mations divided by 10. If it is desired not to use the straight line segment scheme to 
infer the subband levels, it is only necessary to recode as desired the two statements 
defining these variables (and delete the statements defining the line segment slopes). 

Although units are specified in the listing of program BASPAT, the units are, in 
fact, arbitrary except that the input data are required to be in decibels. The units 
specified for temperature, humidity, and distance are dictated by the subroutine 
ATMAT, which is called to obtain discrete-frequency attenuation. 


SUBROUTINE BASPAT (SLI , El ,1 BNF, NF r NB,T, EH, IDK , RO, SLO, OASPL) 

C /BASPAT - BAND SPECTRUM ATTENUATION/ 

C *********************************** 


C * * 

C * ADJUSTS FRACTIONAL-OCTAVE BAND SPECTRUM ANALYSIS SOUND PRESSURE * 

C * LEVEL DATA FOR THE EFFECT OF ATMOSPHERIC ATTENUATION. SINGLE- * 

C * VALUED BAND ATTENUATIONS ARE OBTAINED BY NUMERICAL INTEGRATION * 

C * OF THE DISCRETE-FREQUENCY ATTENUATION EFFECT ON SUBBAND * 

C * RELATIVE LEVELS. SUBBAND RELATIVE LEVELS ARE OBTAINED FROM A * 

C * STRAIGHT LINE REPRESENT ATION OF THE INPUT DATA AS A FUNCTION OF * 

C * THE LOGARITHM OF FREQUENCY. BAND ATTENUATIONS VARY NONLINEARLY * 

C * WITH DISTANCE. IDEAL FILTERS ARE ASSUMED. AIR- AIT ENU ATED * 

C * INPUT DATA YIELD LOSSLESS RESULTS; LOSSLESS INPUT DATA YIELD * 

C * AIR-ATTENUATED RESULTS. * 
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c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* 


♦ 


* SLI(NF) 

* 

* El 

♦ 

* IBNF 

* NF 

* N B 

* T 

* EH 

* IDK 

* 

♦ 

* EO 

* 

* SLO (NF) 

♦ 

* 

* OASPL 

♦ 


INPUT ABE AY OF FRACTIONAL -OCT AYE BAND SOUND PRESSURE * 


LEVELS, DECIBELS * 
PROPAGATION DISTANCE ASSOCIATED WITH INPUT DATA, * 
METERS * 
INTERNATIONAL BAND NUMBER OF FIRST FREQUENCY BAND * 
NUMBER OF FREQUENCY BANDS * 
SIZE OF BANDS, 1/NB-OCTAVE * 
AMBIENT TEMPERATURE, DEGREES KELVIN * 
AMBIENT RELATIVE HUMIDITY, PERCENT * 
INPUT DATA KIND INDICATOR * 
IDK=G DESIGNATES AIR- ATTENUATED INPUT DATA, *' 
OTHERWISE LOSSLESS INPUT DATA ASSUMED * 
PROPAGATION DISTANCE REQUIRED FOR OUTPUT DATA, * 
METERS * 
OUTPUT DATA CORRESPONDING TO INPUT DATA ADJUSTED * 
FOR INVERSE-SQUARE ATTENUATION AND ATMOSPHERIC * 
ATTENUATION * 


OVERALL SOUND PRESSURE LEVEL OF OUTPUT DATA, DECIBELS * 

* 


* CALLS SUBROUTINE ATMAT * 

♦ * 
*************** ************** ****** 


DIMENSION SLI(NF) ,SLOPE{30,2) ,CL(5) ,CR(5) ,SLO(NF) 
M=3 

SIGN=-0. 1 

IF(IDK.EQ.C)SIGN=+0. 1 
R=RO 

IF(IDK.EQ.O) R-RI 
RSQ=2G.0*ALCG10 (RI/RO) 

CB=0. 3/FLOAT (NB) 

CN=2*M+1 
NLS=NF-1 
DO 1 1=1, NLS 

SLOPE (I, 2) = (SLI (1 + 1) -SLI (I) ) /(10.0*CN) 

1 SLOPE (1+1 , 1 ) = SLOPE (1,2) 

SLOPE (1,1) =SLOP E (1,2) 

SLOPE (NF , 2) -SLOPE (NF, 1) 

DU 2 J=1,M 

CL ( J ) =10.0** ( CB/CN* FLOAT (- J) ) 

2 CR (J) =10.0** (CB/CN*FLOAT (J) ) 

0 A SPL=0 • 0 

DO 4 1=1, NF 

FC=1 0.0** (CB* FLOAT (I + IBNF-1 ) ) 

CALL ATMAT (T, RH,R,FC, AC) 

SUMN= 1 0.0** (SIGN* AC) 

SUM D= 1.0 
DO 3 J=1 , M 
F1 = FC*CL (J) 

F2=^FC*CR (J) 

SLI =SLOPE (I, 1) * FLOAT (-J) 

SL2=SLOPE (1,2) *F LGAT ( J) 

CALL ATMAT (T,RH,R, FI ,A1) 

CALL ATMAT(T,RH,R,F2,A2) 

SUMN = SUMN+10.0** (SLI +SIGN* A 1 ) + 1 0 . 0** (SL2+S IGN* A2) 

3 S0MD=SUMD+10. 0**SL1+10. 0**SL2 

SLO (I) =SLI (I) +10.0*ALQG10(SU MN/S U MD) +RSQ 
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JO-O JO-0- 

4 OASPL=OASPL+ (SLO(I)/ 4*3) 

OASPL= 10 . 0* ALOG 1 0 (OASPL) 

RETURN 

END 


Subroutine BAS LA T 


Subroutine BASLAT (BAnd Spectrum Layered ATtenuation) is an alternate version 
of BASPAT that is expanded to allow calculations for a layered atmosphere. The com- 
ments in the text of the preceding section apply to BASLAT as well. 

It should be noted that BASLAT differs only in that it allows the discrete-frequency 
attenuation along the path to be obtained as the sum of incremental attenuations. An- 
other way to accomplish the same result is to use BASPAT as written and to modify the 
program for discrete-frequency attenuation (ATMAT) so that it returns a result ob- 
tained by summation along the path. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* 

* 

* 

* 

* 

* 

♦ 

* 

* 

* 

* 

* 

* 

* 

* 

♦ 

* 

* 

* 

♦ 

* 

* 

♦ 

* 

* 

* 

* 

* 

* 

* 

* 

* 


SUBROUTINE BASLAT (SLI , IDK, IBN F, NF, NB , NL,T,fiH, RI NC ,R L, S LO ,OAS PL) 
/BASLAT - BAND SPECTRUM LAYERED ATTENUATION/ 

********************************* 


THIS IS AN ALTERNATE VERSION OF BASPAT THAT HANDLES INCREMENT- * 
ALLY VARIABLE ATMOSPHERIC PROPERTIES ALONG THE PROPAGATION * 

PATH, WHICH MIGHT BE REQUIRED BY A LAYERED ATMOSPHERE MODEL, * 

FOR EXAMPLE. IF A SINGLE INCREMENT (LAYER) IS SPECIFIED * 

(I.E., NL=1 ) , THIS PROGRAM PERFORMS EXACTLY AS BASPAT. * 

* 

SLI (NF) INPUT ARRAY OF FR ACTIONAL -OCT AVE BAND SOUND PRESSURE * 

LEVELS, DECIBELS * 

IDK INPUT DATA KIND INDICATOR * 

IDK =0 DESIGNATES AIR- A TTEN UA TED INPUT DATA, * 

OTHERWISE LOSSLESS INPUT DATA ASSUMED * 

IBN F INTERNATIONAL BAND NUMBER OF FIRST FREQUENCY BAND * 

NF NUMBER OF FREQUENCY BANDS * 

NB SIZE OF BANDS, 1/NB-OCTAVE * 

NL NUMBER OF PATH LENGTH INCREMENTS; FOR A HOMOGENEOUS * 

ATMOSPHERE, SET NL=1 * 

T (NL) AIR TEMPERATURE IN EACH PATH INCREMENT, DEGREES * 

KELVIN * 

RH (NL) RELATIVE HUMIDITY IN EACH PATH INCREMENT, PERCENT * 

RINC(NL) PROPAGATION PATH LENGTH INCREMENTS FOR THE AIR- * 

ATTENUATED DATA, METERS * 

RL PROPAGATION DISTANCE IN METERS FOE THE LOSSLESS DATA * 

SLO (NF) OUTPUT DATA CORRESPONDING TO INPUT DATA ADJUSTED * 

FOR INVERSE-SQUARE ATTENUATION AND ATMOSPHERIC * 

ATTENUATION * 

OASPL OVERALL SOUND PRESSURE LEVEL OF OUTPUT DATA, DECIBELS * 

* 

CALLS SUBROUTINE ATMAT * 

* 


****************** **************** 

DIMENSION SLI (NF) ,SLOPE(30,2) ,CL(5) ,CK(5) ,SLO(NF) 

DIMENSION T (NL) , RH (NL) , R INC (N L) 

M=3 
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SIGN=-0. 1 

IF (IDK.EQ.O) S IGN =+0 . 1 

RA=0.0 

DO 1 K=1,NL 

1 RA=RA+RINC (K) 

RSQ= (10.0*SIGN) *20.0*ALOG10 (RA/RL) 

CB = 0 . 3/F LOAT(NB) 

CN=2*M+1 
NLS=NF-1 
DO 2 1=1, NLS 

SLOPE (I,2) = (S LI (1+1) -SLI(I) ) /(I 0.0* ON) 

2 SLOPE (1+1 , 1) =SLOPE (1,2) 

SLOPE (1, 1) =SLOPE (1 ,2) 

SLOPE (NF , 2) =S LOPE (NF , 1) 

DO 3 J=1,M 

CL ( J ) =10.0** (CB/CN* FLOAT (-J) ) 

3 CR (J) =1 0.0** (CB/CN*FLOAT(J) ) 

OASPL=0 . 0 

DO 7 1=1, NF 

FC=1 0.0** (CB* FLOAT (I + IBNF-1 ) ) 

A C= 0 . 0 
DO 4 K= 1 , N L 

CALL ATMAT (T (K) , RH ( K) ,RINC(K) , FC , DELTA) 

4 AC=AC +DE LTA 
SUMN=10.0** (S IGN* AC) 

S UMD= 1 . 0 

DO 6 J= 1 , M 
F1 = FC*CL (J) 

F2= FC*CR (J) 

SL1= SLOPE (1,1) * FLOAT (-J) 

SL2 = SLOPE (I, 2) * FLOAT (J) 

A 1 = 0 . 0 
A2 = 0 .0 
DO 5 K= 1 , NL 

CALL ATMAT (T ( K) ,RH (K) ,RINC(K) ,Fl,DiiLTA1) 

CALL ATMAT (T ( K) , RH (K) , RINC(K) , F2 , DELTA2) 

A 1 = A 1 +DE LT A 1 

5 A2= A2 +DLLT A2 

SUMN = SUMN+10.0** (SLl +SIGN*A1) + 1 0 . 0 * * ( S L 2+S IGN * A2) 

6 SUMD=SUMD+1Q. 0**SL1 +10. 0**SL2 

SLO (I) =SLI (I) ^ji.0*ALOGl0 (SUMN/SUMD) + RS Q 

7 OASPL = OASPL+ ^5** (SLO (I)/ 7<dflO*0 
OASPL= 10. 0*ALOG 10 (OASPL) 

RETURN 

END 


Subroutine ATMAT 

Subroutine ATMAT (for ATmospheric ATtenuation) is called by BASPAT and 
BASLAT to obtain total-path discrete-frequency atmospheric attenuation for a homoge- 
neous atmosphere. The substance of this program is taken from reference 2. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* 

* 

* 

* 

* 

* 

* 

♦ 

♦ 

♦ 

♦ 

* 

♦ 

* 

* 

* 

* 


SUBROUTINE ATfiAT (T , RH , DIST ,FREQ, ATT) 

/ATMAT - ATMOSPHERIC ATTENUATION/ 

********************************** 

* 

COMPUTES ATMOSPHERIC ATTENUATION IN DECIBELS FOR SPECIFIED 
TEMPERATURE, RELATIVE HUMIDITY* DISTANCE* AND FREQUENCY 
FOLLOWING BASS AND SHIELDS* NASA CR-2760, 1977. 

APPLICABLE TO FREQUENCIES FROM 50 TO 100000 HERTZ* TEMPERATURES 
FROM 255.4 TO 310.9 DEGREES KELVIN* RELATIVE HUMIDITIES FROM 
ZERO TO ONE HUNDRED PERCENT. 


********** 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

*********************** 


T TEMPERATURE, DEGREES KELVIN 

RH RELATIVE HUMIDITY, PERCENT 

DIST DISTANCE, METERS 

FREQ FREQUENCY, HERTZ 

ATT ATTENUATION, DECIBELS 

P AMBIENT PRESSURE IN ATMOSPHERES, MAY BE MADE A 

CALLING VARIABLE (SEE CR-2760) 


P = 1 .0 
Tl=T/293. 

T01 =273 . 16 

PS =10. 79586* (1.0- TO 1/T) - 5. 028 08*ALOG 1 0 { T/T0 1 ) 

1 + 1. 50474E-4* (1.0-10.0** (-8 . 29 692* ( (T/TO 1) -1.0) ) ) 

2 +0.42 87 3 E-3* (10.0** (4.76955* (1.0- (T0 1/1) ) ) - 1 . 0) -2 . 2 1 95 983 
PS=1G.0**PS 

H=RH* PS/P 

FRO=P* (24. 0+4 ,41E+4*H*(0.05+H)/(0.391 + H) ) 

FRN=P/SQRT (T1 ) * (9 . 0 +350 . 0* H* EXP ( -6 . 14 2* { (1.0/T1) **0.333-1.0) ) ) 
ALPHA=SQRT ( T 1 ) * FREQ** 2/P* ( 1 . 84E- 11 + 2.19 13E- 4/Tl*P* (223 9 . 1/T) **2 
1 *EXP (-223 9 . 1/T) / (FRO+ (FREQ* *2/FRO) ) +8. 16 19E-4/T1*P* (3352. 0/T) **2 - 

2* EXP (-3 352. 0/T) / (FRN+ (FREQ**2/FR N) ) ) *8.686 
ATT=ALPHA* DIST 
RETURN 
END 
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TABLE I. - COMPARISON OF NUMERICAL INTEGRATION RESULTS 


FROM PRESENT METHOD WITH THOSE OF REFERENCE 2 


Interna- 

tional 

band 

number 

1/3-Octave 
band nomi- 
nal center 
frequency, 
Hz 

1/3 -Octave band sound pressure levels at 
20 meters from source 

Original 

lossless 

Air -attenuated at 
293 K, 50 percent 
relative humidity 

Inverse calculation 
to lossless 

Present 

method 

Reference 2 

Present 

method 

Reference 2 

36 

4 000 

40.0 

39.5 


40.0 


37 

5 000 

37. 0 

36.2 

36.2 

37.0 


38 

6 300 

34.0 

32.8 

32. 8 

34.0 

34 

39 

8 000 

31.0 

29. 1 

29.2 

31.0 

31 

40 

10 000 

28.0 

25. 1 

25.2 

28.0 

28 

41 

12 500 

25.0 

20.6 

20.8 

25.0 

25 

42 

16 000 

22.0 

15.4 

15. 5 

22.0 

22 

43 

| 20 000 

19.0 

9.3 

9. 5 

19.0 

19 

44 

25 000 

16. 0 

1.9 

2.2 

16.0 

16 

45 

31 500 

20. 0 

. 1 

.4 

20.0 

20 

46 

40 000 

24.0 

-2.5 

-2.4 

24.0 

24 

47 

50 000 

28.0 

-6. 1 

-5.6 

28.0 

28 

48 

63 000 

32. 0 

-10.9 

-10. 2 

32.0 

32 

49 

80 000 

36.0 

-17. 5 

-17. 1 

36.0 

-- 

50 

100 000 

40. 0 

-27.3 


41.0 

-- 


i^ 1 Band, n ° 2m + 1 subbands, -m ^ j < + m 


Subband 



Log (frequency) — ► 


Figure 1. - Delineation of fractional-octave bands and subbands on logarithmic fre- 
quency scale. 







f 


t 



Log (frequency) — ► 


Figure 3. - Scheme to infer subband level distribution across each band of fractional -octave 
band analysis using straight-line half-segments. Solid symbols denote band levels; open 
symbols denote subband levels; m * 2, n - 2m + 1 - 5, for this example. 



Figure 2. - Relative slope and amplitude relationships between band levels and power spectral density, 


Error in lossless data, [a (f c ) • r 


i th Band 



Log (frequency) — ► 


Figure 4. - Subband level distribution relative to subband 
level at band center frequency. 



1/3-Octave band center frequency, kHz 

(a) Frequency range, 2 to 20 kilohertz for 500-meter (b) Frequency range, 10 to 100 kilohertz for 10-meter 

propagation distance. propagation distance. 

Figure 5. - Error in lossless 1/ 3-octave band data obtained from air-attenuated data at 298. 15 K, 70 percent relative humidity 
by using center frequency attenuation a{f c ) - r instead of integrated results from present method Aj(r). 
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Error in lossless data, A:{r 



1/ 3-Octave band center frequency, kHz 

(a) Frequency range, 2 to 20 kilohertz for 500-meter (b) Frequency range, 10 to 100 kilohertz for 10-meter 

propagation distance. propagation distance. 

Figure 6. - Error in lossless 1/3-octave band data obtained from air -attenuated data at 298. 15 K, 70 percent relative humidity 
by assuming band attenuation is proportional to distance. Attenuation for each unit distance obtained by numerical in- 
tegration and proportioned to actual propagation distance to obtain approximate results A j(r u ) - r/r u . Resulting differences 
from integrated results for actual distance Aj(r) are shown. 
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